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Abstract 


Hingeless rotor frequency response calcula- 
tions are obtained by applying a generalized har- 
monic balance to the elastic blade flapping equa- 
tions. Nonuniform, unsteady Induced flow effects 
are Included by assuming a simple three-degree- 
of -freedom description of the rotor wake. Results 
obtained by using various models of elastic blade 
bending and induced flow are compared with exper- 
imental data obtained from a 7.5-ft diameter wind 
tunnel model at advance ratios from 0.0 to 0.6. 

It is shown that the blade elasticity and nonuni- 
form, unsteady induced flow can have a signifi- 
cant effect on the transient response character- 
istics of rotor systems. Good correlation be- 
tween theory and experiment Is obtained by using: 
(I) a single rotating mode shape description of 
the elastic blade bending, (li) an empirical form- 
ula for the quasi-steady induced flow behavior, 
and (ill) the apparent mass terms from potential 
flow for the unsteady induced flow characteris- 
tics. 
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two-dimensional lift-curve slope, 
rad“l 

harmonics of jth flapping mode 
number of blades 
tip loss factor 
blade chord, ft 

steady value of thrust coefficient, 
steady thrust /pirS^R 1 * 
harmonic perturbation of thrust 
coefficient 

harmonic perturbation of roll moment 
coefficient * roll moment /p*Q 2 R 5 , 
positive advancing blade down 
harmonic perturbation of pitch 
moment coefficient * pitch moment/ 
pirfl 2 R 5 , positive nose up 
dimensionless flapping hinge offset 
dimensionless radius of pocket 
cutout 

rotor blade bending stiffness, 
lb-ft 2 

generalized response vector 
aerodynamic and inertial forces per 
unit blade span, lb /ft 
nondlmensional harmonics of inertial 
forcing function, Eq. (12) 
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number of flap bending modes 
nondimenslonal apparent mass and in- 
ertia of impermeable disk 
control feedback matrix 
nonuniform induced flow matrix 
empirical value for quasi-steady 
portion of [L] 

rotor blade mass distribution, slug/ft 
nondimenslonal blade parameters 
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rotor response matrix open loop, 
closed loop 

elemental apparent mass, slugs 
elemental mass flow, slugs /sec 
apparent mass of air, slugs 
index referring to harmonic number 
number of azimuthal harmonics 
null matrix 

first flap frequency divided by H 

generalized coordinates 

steady values of qj 

rotor blade radius coordinate, ft 

rotor blade radius, ft 

blade root moment, ft-lb 

blade parameter 
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generalized control vector 
< ®o®s®c®oSs8c^o^s^c > 

perpendicular, tangential components of 
air speed in undeformed blade coordi- 
nate system, ft /sec 

freestream airspeed perpendicular and 
parallel to rotor shaft (V® positive 
down), ft/sec 

induced flow parameter * 

[U 2 + X(\ + v)]/(y 2 +1 2 ) 1/2 

blade root shear, lb 

rotor blade flap deflection, ft 

frequency transform, Eq. (23) 

physical control vector <0 o 0 s 0 c z^a> 

control coupling matrix 

hub plunge deflection divided by R, 

positive down 

hub pitch angle, positive nose up, 
rad 
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Lock number, 1/myy 

equivalent Lock number, Eq. (35) 

blade pitch angle - 5 + (0 O + 0 S 

sin ^ + 0 C cos ij/)e i{ 4 

steady collective pitch angle 

rotor pitch perturbations 

total inflow (including induced flow) 

= X + ( x 0 + x s | sin i, 

+ X c ^ cos e iu ^ 

steady inflow ratio ■ V ro /ftR + v 
inflow perturbations (including in- 
duced flow), Eq. (10) 
advance ratio « Uw/ftR 
total induced flow * 

v + ( V D + v 8 | sin ill 

, r \ iijjiJj 

+ v c — cos ij; j e 

induced flow due to steady rotor 
thrust 

induced flow perturbations 

air density, slug/ft 3 

rotor solidity, bc/rR 

induced flow time constants, rad”l 

hub roll angle, positive advancing 

blade down, rad 

orthogonal functions 

rotor blade azimuth position, non- 

dimensional time, rad 

excitation frequency divided by ft 

rotor blade angular velocity, 

rad /sec 

3/3r 

3/a^ 


The dynamic response characteristics of 
hingeless rotors are dependent upon the distrib- 
uted structural properties of the rotor blades, 
the local aerodynamic properties of the blade 
sections, and the detailed description of the 
aerodynamic environment. It is generally be- 
lieved, however, that reasonable predictions of 
rotor thrust and moments at low lift can be ob- 
tained by using some appropriately simplified 
models for the blade structure, section aero- 
dynamics, and inflow distribution. The develop- 
ment of these simplified rotor models is useful 
for gaining Insight into the basic dynamic mech- 
anisms of rotor response. Detailed calculations 
of dynamic airloads, necessary for many applica- 
tions, are usually too complex for use in basic 
dynamic research or preliminary design calcula- 
tions. 


The formulation of a minimum complexity 
rotor response model is the subject of several 
recent papers. One area of interest is the ef- 
fect of mode shape and mode number on rotor 
flapping response. Shupe 1 addresses the effects 
of the second flap mode, Ormlston and Peters 2 
compare various mode shape models for first and 
second flap modes, and Hohenemser and Yln^ con- 
sider the effect of using rotating rather than 
nonrotating modes as generalized degrees of 
freedom. The fundamental conclusion, as 


clarified in Reference 3, is that for y < 0.8 a 
single rotating mode shape is adequate for model- 
ing the steady rotor flapping response. 

A second area of interest is the effect of 
induced flow perturbations on rotor flapping 
response. In Reference 1, a simple momentum 
theory predicts a significant effect of induced 
flow on steady rotor response. In Reference 2, a 
comparison of steady experimental and theoretical 
results indicates that, although there is a sig- 
nificant effect due to induced flow, momentum 
theory is inadequate for predicting this effect in 
forward flight. Alternate induced flow models are 
introduced and compared with the data, but no 
clear choice for the best model is found. In 
Reference 4, an unsteady momentum theory is used 
in hover to improve correlations with experimental 
frequency response data. 

The work in References 1 through 4 indicates 
that a minimum complexity analytic model for 
rotor dynamics must include appropriate degrees of 
freedom for both structural and induced flow per- 
turbations (certain flight dynamics programs 
presently include a simplified dynamic treatment 
of the induced flow^). Unfortunately, while some 
success has been achieved using simple models of 
the rotor induced flow in hover, a completely 
satisfactory induced flow model for forward flight 
has not been found, not even for the condition of 
steady response. In addition, neither the phy- 
sical values of the induced flow time constants 
nor the frequency range in which they are impor- 
tant is known. The unsteady behavior of the in- 
duced flow contributes directly to the low fre- 
quency rotor control characteristics and to the 
coupled rotor /fuselage aeroelastic stability. In 
particular, induced flow perturbations contribute 
to the rotor damping available in pitch and roll 
(which is important for ground and air resonance 
calculations) . It is consequently important to 
understand the dynamic characteristics of the in- 
duced flow. 

The purpose of this paper is to provide ad- 
ditional insight into the question of rotor 
structural and induced flow modeling. To this 
end, experimental rotor frequency response data 
in hover and in forward flight are compared with 
theoretical results that are calculated by using 
several different models for the elastic blade 
bending and Induced flow. The frequency response 
data provide a broad base of comparison so that 
the effects of mode shape and Induced flow model 
can be clearly determined throughout the fre- 
quency range of interest. 

Basic Equations 


Analysis 

The mathematical technique used here is a 
further generalization of the harmonic balance 
approach of Reference 2. In addition to an arbi- 
trary number of bending modes (with an arbitrary 
number of azimuthal harmonics for each mode) , the 
generalized harmonic balance allows for a 
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rational treatment of reversed flow aerodynamics 
and the possibility of harmonically oscillating 
control inputs. 


The linear equation of motion for the de- 
flection of an elastic beam subject to distrib- 
uted aerodynamic and inertial loadings F(r,4>) 
and G(r,*) is® 

(Elm”) + + fl 2 ^mrw f - w" J mr drj 


are obtained for the a* and bj n . Solution of 
these equations, followed by a substitution of Eq. 
(5) into Eqs. (2) and (3), results in the phase 
and magnitude of all desired harmonics of the 
flapping deflections and hub forces and moments. 

Blade Loading 

The aerodynamic loading of each blade is 
given by 


• F(r,*) + G(r,^) (1) 

The associated expressions for bending moment and 
shear at the blade root are 

r R 

S(0,*) - J (F + G - mQ^w - mfl^r dr (2) 
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Up 3 Qw + ftRA + ftRpw' cos $ 


( 8 ) 


The blade root bending moment is transformed into 
a stationary coordinate system to yield the pitch 
and roll moment of the rotor. The solution of 
Eq- (1) yields directly the blade deflections, 
and substitution into Eqs. (2) and (3) then 
yields the forces and moments. 

Application of the harmonic balance involves, 
first of all, an orthogonal expansion of w: 

J 

f - % qjMtjCr) (4) 

J-l 


Eq. (8) contains the primary contributions of mode 
shape and induced flow to the flapping equations. 
The details of blade mode shape become important 
as p increases because Up depends upon both 
the blade deflection w and its first derivative 
v f . The induced flow is important because first 
order perturbations to the inflow X create first 
order changes to Up and F. 

Although the inflow is in general a compli- 
cated function of radius and azimuth, as a first 
approximation, the total inflow can be represented 
by 


For the present analysis, the are taken to 

be the exact mode shapes of the rotating beam 
without aerodynamics. Galer kin's method is then 
used to transform Eq. (1) into J ordinary dif- 
ferential equations (with periodic coefficients) 
for the modal coordinates qj.6 When the forcing 
terms contain a steady portion superposed onto 
periodic functions that are modulated by an ex- 
citation frequency oj (cycles per revolution) , 
Floquet's theorem Implies that the qj have a 
solution of the form' 


q j * q j + 
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+ X S R Sin * + X C f COS 


The steady portion of the total inflow X con- 
tains contributions from the freestream velocity 
V m /QR and from the steady induced flow due to 
rotor thrust v. The unsteady inflow components 
X 0 ,X S ,X C contain contributions from harmonic 
plunging ze^, rolling ^e^, and pitching 
ae iaj ^ of the shaft, as well as contributions 
from the unsteady induced flow components 
v 0 ,v s ,v c due to perturbations in rotor thrust and 
moments : 


X * -iwz + v - pa 
o o I 


X g * -iio<j> + v g V 


(10 a-c) 


where q are the steady coning displacements 
and the J aj n and bj n are complex quantities 
indicating the magnitude and phase shift of each 
modulated harmonic of the perturbation response. 
The harmonic balance approach entails substitut- 
ing Eq. (5) into the J ordinary differential 
equations for qj and Retting coefficients of 
like harmonics equal. When n is truncated at 
the highest harmonic of interest N, then 
(2 • N + 1) • J linear algebraic equations 


X = -iwa + v 9 

c c J 

The blade pitch angle 0 is given by 
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where 0 is the steady value of 0 and 
®o* e s» e c are contr °l system perturbations. The 
inflow perturbations X Q ,X S ,X C are assumed to be 
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small compared with unity. This implies that the 
induced flow perturbations v 0 ,v g ,v c and the con- 
trol perturbations 6 0 > e s*9 c *z»^,ot are also small 
quantities yielding linear perturbation equa- 
tions. 

The inertial loading of each blade is given 
by 


moments influence the induced flow. The induced 
flow, therefore, is a feedback loop of Eq. (15), 
causing the uj to depend upon the f j . 

From the standpoint of calculation, it is con- 
venient to express the coupling relation (between 
the generalized controls, the physical controls, 
and the rotor response) in matrix form: 


G = -mfi 2 R^g 0 + g s ^ sin i|i + g c ^ cos i|»J e 1 ^ (12) 

where 

g Q =* U) 2 Z I 

g g * w 2 <f> + 2io)a ) (13 a-c) 

g * -2iw<J> + w 2 a ! 

c / 

The inertial loading is a result of centrifugal, 
Coriolis, and gyroscopic forces which occur in 
the rotating reference frame of the blade due to 
hub motions z,<f>,a in the inertial reference 
frame. 

When Eqs. (6) through (12) are combined and 
appropriately integrated in Eqs. (1), (2), and 
(3), the steady deflections and forces q.i ,5p/aa 
are obtained as linear functions of the steady 
inputs 6, A; and the perturbation blade deflec- 
tions and hub forces and moments are obtained as 
linear combinations of the generalized control 
variables 

<u>=<^eeegggAxx \ (14) 

N O 8 C 6 0 6 S 6 C o 8 C ' K ' 

Although g 0 ,g s ,g c are simply related to the 
shaft motion through Eq. (13), they are retained 
as generalized controls so that the generalized 
controls can be separated into physical, in- 
ertial, and aerodynamic groupings. This will 
facilitate the calculation of rotor response when 
induced flow is included later. 

Interpretation of Results 

The results of the harmonic balance can be 
expressed in matrix form as 

{f} - [M] {u} (15) 

where {f} represents the perturbation harmon- 
ics of thrust, moments, and generalized coordi- 
nates. The elements of [M] , therefore, have 
direct physical significance. They are the par- 
tial derivatives of each of the response harmon- 
ics taken with respect to each of the generalized 
controls u j . The generalized control variables 
are in turn functions of the physical controls 
xj. 


{u} = [Y] {x} + [K] { f } (17) 

Eq. (17) is simply a set of linear equations de- 
scribing: (i) the generalized control perturba- 

tions due to application of the physical controls 
[Y] and (ii) the generalized control perturba- 
tions due to the effect that rotor response has 
on the induced flow [K] . The matrices [Y] and 
[K] will be obtained later by using an appropri- 
ate induced flow model. It follows that the par- 
tial derivatives of the f j with respect to the 
physical controls x^ can be found (including in- 
duced flow effects) from Eqs. (15) and (17). The 
derivative matrix is designated [M T ] and has the 
properties 


(f> * [M']{x} 


(18) 


[M* ] - |I] - CM] £K] 


-1 


[M][Y] 


(19) 


Although the higher harmonics are often necessary 
in the harmonic balance calculation of [M] , the 
subsequent calculation of [M f ] by Eq. (19) may 
be performed for only those response and inflow 
harmonics of interest. In this paper, five har- 
monics are used in the calculation of [M] , but 
only first harmonics are retained in Eqs. (18) and 
(19) , so that the f j are taken to be 
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Induced Flow 


Form of Induced Flow Model 


A useful form of the induced flow model is 
given by ^ 
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aerodynamic only 


( 21 ) 


Although not completely general, Eq. (21) can 
accommodate a variety of induced flow models. 
Only aerodynamic contributions are included on 
the right-hand side, because they are the only 
loads which produce reaction forces on the rotor 
wake. Using Eqs. (2) and (3), these aerodynamic 
forces and moments can be expressed in matrix 
form as 


<X> *< 9 O 6 8 e c Zl(,0l > <16) 

as evidenced in Eqs. (10) and (13). 

The generalized control variables ua are 
also coupled to the f j , because the thrust and 
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With the induced flow v described by Eqs. (21) 
and (22), the inflow relation follows directly 
from Eqs. (10), (13), and (21). The matrices 


m 

and [K] 

of Eq. (17) may then be 

as 
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Eq. (24) represents the control coupling be- 
tween the physical controls xj and the general- 
ized controls u j . The presence of [L] in this 
matrix Indicates that the X's are indirectly 
coupled (through the induced flow), as well as 
geometrically coupled [Eqs. (10) and (13)] to the 
rotor plunge, pitch, and roll motions. Eq. (25) 
represents the induced flow caused by the de- 
pendence of X upon the thrust and first harmon- 
ic flapping. If a suitable approximation to the 
inflow can be modeled in the form of Eq. (21), 
then Eqs. (24) and (25) may be substituted 


directly in Eq. (19) to obtain the complete rotor 
response to physical control inputs. 

Unsteady Momentum Theory 

An approximation of the induced flow that is 
suitable for Eq. (21) can be obtained as an ex- 
tension of the momentum theory used in Reference 
2. The differential force on an elemental area of 
rotor disk is written as 

dF « 2KRvdm + n 2 Rvdm (26) 

where 2URv is the total change in velocity nor- 
mal to the disk, dm is the differential mass flow 
through the element, dm is the v apparent mass 
associated with the flow, and v is the time 
derivative of v in the nonrotating system. The 
differential mass flow relation 


dm * pQR Vp 2 + X 2 r dr d^ (27) 

can be used to integrate the first term of Eq. 

(26) over the disk to obtain a quasi-steady in- 
duced flow relation for rotors that have combined 
conditions of thrust and forward speed. The eval- 
uation of the second term In Eq. (26) (the un- 
steady effect) requires the additional knowledge 
of the apparent mass dm associated with the flow. 

An approximation to the apparent mass terms 
of a lifting rotor can be made in terms of the 
reaction forces (or moments) on an impermeable 
disk which is instantaneously accelerated (or ro- 
tated) in still air. This approximation was used 
in Reference 8, giving good agreement with trans- 
ient thrust measurements for an articulated rotor. 
The reactions on such an impermeable disk are 
given from potential flow theory in terms of el- 
liptic Integrals which are evaluated in the lit- 
erature. ^ They result in apparent mass and in- 
ertia values 


”A “ f PR3 ' *A * If PR5 < 28 > 

(For v * v 0 r/R, a radial velocity distribution, 
m A * pR 3 .) These values represent 64 percent of 
the mass and 57 percent of the rotary inertia of 
a sphere of air having radius R. It is empha- 
sized that they are only approximations to the 
actual values for a lifting rotor. 

Using this approximation, the steady ind uced 
flow equation and the unsteady induced flow per- 
turbation equations can be derived from Eqs. (26) 
through (28) : 

2v Vy 2 + X 2 = Qj, 

K v + 2vv * C_ 
mo o T 

Vs + I VV S - "S. 

Vc + I vv c - -Si 


(29a) 

(29b-d) 
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where 


V= Jti + X < X + il <29e) 

\/u 2 + * 2 

and 


m A fl 

K = — = — = 0.8488 ) 

m pirR 3 3it / 

> (30a-b) 

K 3 — = — = 0.1132 ' 
pirR 5 45tt 


Eq. (29a) expresses the nonlinear relation be- 
tween the steady thrust and the steady induced 
flow v. Eqs. (29b-d) are then the linear per- 
turbation equations for small changes in thrust, 
moments, and induced flow. In order for the per- 
turbation equations to be valid, it is assumed 
that v 0 > v 8 ,v c are much smaller than (u 2 + A 2 )^. 


The time constants associated with the induc- 
ed flow model in Eq. (29) are 


m 

2 v 

= 0.4244/v 

( f or v o ) i 


2K i 

v 

= 0.2264 /v 

(f° r v s ,v c ) ) 

(31a-b) 


In Reference 4, the steady induced flow v and 
the time constant for v s ,v c are obtained by 
correlating experimental hover frequency response 
data. Two operating conditions are considered, 
and the best fit in these cases is found to be 
v = .014, t s - 8 (with 0 = 2 °) and v = .028, 

Tg ** 4 (with 6 * 8 °). From the v values indi- 
cated for these cases, it can be shown that each 
t s implies the same value of Kj « 0.112. Thus, 
there is some experimental evidence that the po- 
tential flow value K x = 0.113 is approximately 
valid. 


By assuming simple harmonic motion, Eqs. 
(29b-d) can be brought into the form of Eq. (21), 
yielding the components of [L] for unsteady 
momentum theory. 


qa 

2 v + K iui 
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[L] 


v /2 + K^irn 


(32) 


,/2 + K t 


io) 


(L 22 and L 33 differ by a factor of 4/3 from Ref- 
erence 2 , because v s and v Q are taken uniform 
with r in that reference, whereas they are 
taken linear with r here.) The matrix [L] 
from Eq. (32) may now be substituted into Eqs. 


(24), (25), and (19) to obtain the rotor response 
that includes inflow. 


Empirical Model 


Experimental data have shown that momentum 
theory, although particularly simple to use, is 
qualitatively inaccurate for certain steady re- 
sponse derivatives in forward flight. Reference 
2 introduces an alternate induced flow model for 
forward flight in which the elements of [L] 

(with to * 0 ) are chosen to give the best fit of 
experimental response data for several configura- 
tions at conditions of near zero lift. If this 
empirical inflow model, [L E ] * is taken for the 
quasi-steady portion of the induced flow law, and 
if the theoretical apparent mass terms (from po- 
tential flow) are taken as a model for the un- 
steady portion of the induced flow law, then a 
complete induced flow equation can be expressed as 



[C T /aa 
Jc^/oa^ 
Jyaa! 


(33) 


The assumption that the apparent mass terms 
may be superposed on the quasi-steady terms is not 
rigorous, but it can be considered analogous to 
unsteady wing theory in which the apparent mass 
terms are theoretically independent of the free- 
stream velocity. Under the superposition assump- 
tion, the empirical inflow model modified for the 
unsteady case is 
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0 
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~ K I 
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aa 

l e 

_0 
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- K I- 
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(34) 


Although this particular formulation of [L] is 
valuable for predicting the effects of induced 
flow, ultimately a more consistent formulation of 
[L] should be made, as discussed in Reference 10. 


Equivalent Lock Number 

Another method of accounting for the unsteady 
induced flow is the use of an equivalent Lock num- 
ber y* , which can be derived from a single har- 
monic balance of the root moment equation: 



1 _ 

1 + 8 v/aa + 16K I iu)/ aa 


(35) 


Although this approach is not a completely con- 
sistent treatment of the induced flow, since it 
does not give an exact harmonic balance of the 
blade flapping and thrust equations, it yields 
results which are nearly the^same as those ob- 
tained from momentum theory. 
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The practical use of Eq. (35) is somewhat 
limited because of the inaccuracies of momentum 
theory in forward flight, but a y* approach is 
nevertheless a valuable conceptual tool for under- 
standing the effects of induced flow. In particu- 
lar, Eq. (35) shews that one effect of induced 
flow perturbations is to decrease the effective 
Lock number (i.e., decrease the aerodynamic ef- 
fectiveness). This decrease is most pronounced at 
low values of v (i.e., low y and 9 ) and low 
values of u>. For example, rotor roll moment is 
plotted in Figure 1 for two values of 9 and com- 
pared with the value from elementary theory 
(steady induced flow only, induced flow perturba- 
tions neglected, equivalent to 11m 9 -► 00 ). The 
curves for 0 * 0, 0.05 result in values of roll 
moment well below the elementary value. 



e s , rod 


Figure 1. Effect of induced flow on steady rotor 
response in hover, y * 0, u> * 0, 
o - 0.1, a - 2x, p - 

The effect of induced flow is most pro- 
nounced in the response derivative (the slope of 
the response curve at 6 S - 0). For p - ®, the 
derivative is given by 


3 [c L / a.J/a e s l =-&f(l + 3/2u*) (36) 

' 8 

indicating that y*/y < 1 results in a reduction 
of the roll moment response (or control power) 
from the elementary value. When the rotor Is in 
hover with no lift (v * 0), a quasi-steady per- 
turbation of 6 8 (u) • 0) results in no response 

because of the aero slope of the curve in Figure 
1. The mathematical justification for the van- 
ishing response derivative can be seen in Eqs. 

(35) and (36) . With u - v - 0, y*/y and the 
response derivative must equal zero. As ? in- 
creases, however, v and v increase so that 
Y*/Y approaches unity and the derivative ap- 
proaches -1/16, as illustrated in Figure 2. With- 
in the practical range of thrust coefficients, 
however, the response derivative never recovers 
more than about 80 percent of the elementary 
value. Eq. (35) also Implies that increasing ad- 
vance ratio (which increases v) will result in a 
partial recovery of y*/y (and of the response 
derivative). This recovery is evident in 
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Figure 2. Effect of induced flow on steady rotor 
response derivatives in hover, y * 0, 
w = 0, a = 0.1, a ® 2 tt, p - 

Figure 3, where the roll response is given versus 
y; but no more than 90 percent of the elementary 
value is reached in the practical range of thrust 
and advance ratio. 



Figure 3. Effect of induced flow on steady rotor 
response derivatives in forward 
flight, id ■ 0, o » 0.1, a - 2ir, p - 


The unsteady terms (apparent inertia K-) 
also bring y*/y closer to unity, as seen by the 
role of K x in Eq. (35). This recovery with fre- 
quency is illustrated in Figure 4, where, as u 
becomes large, the response derivative approaches 
the elementary value of -1/16. The rate at 
which the response approaches -1/16 is dependent 
upon the magnitude of the apparent Inertia K x . 
Large values of K x result in a rapid return to 
the elementary value, and small values of K x re- 
sult in a slow return. For K x ■ 0.1132 and 
u < 0.3, the unsteady terms provide only small 
contributions to the response. Thus, the quasl- 
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Figure 4. Effect of induced flow time constant on 
rotor frequency response derivatives, 
y * 0, a = 0.1, a = 2 it , p = », 
v * X = 0.05. 


steady theory (with K x = 0) would be adequate in 
this range. In the frequency range 0.3 < go < 

1.0, the unsteady terms have a more significant 
effect. Above co = 1.2, the total effect of in- 
duced flow diminishes so that the elementary 
theory and the unsteady theory give similar re- 
sults; but the quasi-steady theory (with K x - 0) 
is in considerable error in this region. 

The frequency range in which unsteady in- 
duced flow is important_is also dependent upon the 
thrust or mean inflow v as shown in Figure 5. 

For low values of v, the unsteady effects domi- 
nate at low frequencies; and for large values of 
v, the unsteady effects are delayed into the 
higher frequencies. This effect is implicit in 
Eq. (35) and is a direct result of the inverse 
dependence of_time constant upon v, Eq. (31). 
Thus, a low v implies a slow induced flow re- 
sponse; and a high v implies a rapid induced 
flow response. Equation (35) shows that advance 
ratio (which also increases v) has a similar ef- 
fect on the induced flow behavior. It follows 
that the relative importance of the unsteady and 
quasi-steady nonuniform induced flow terms de- 
pends upon both the rotor operating conditions 
and the frequency range of interest. 



UJ 


Figure 5. Effect of unsteady induced flow on 

rotor frequency response derivatives, 
y ® 0, a = 0.1, a = 2tt, p = °°, 

K x = 0.1132. 
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In Figure 6, the relative importance of 
these terms is presented qualitatively through a 
chart of the operating regimes in which (for no 
induced flow or quasi-steady induced flow) |y*| 
differs by less than 10 percent from the unsteady 
value. This is a subjective criterion and is 
merely intended to illustrate the trends with 
thrust, advance ratio, and frequency. Four 
regions are defined: (i) at high go and v, 

induced flow effects are small and either the 
elementary or quasi-steady approximation is ade- 
quate; (ii) at high go and low v, although in- 
duced flow effects are small (no induced flow 
being a good approximation) , the quasi-steady 


Figure 6. Regions of validity for steady (no 
induced flow perturbations) quasi- 
steady (K x = ~ 0) , and unsteady 

(K x = 0.1132, Kjh - 0.8488) induced 
flow models based on Eq* (35), 

a - 0.1, a = 2 it . 

theory alone will be in error; (iii) at low go 
and high v, the opposite is true (i.e., the 
quasi-steady nonuniform theory is required, where- 
as neglecting induced flow results in error) ; and 
(iv) for low go and v, complete unsteady theory 
is required. 
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Comparison of Theory and Experiment 

The experimental data used in the following 
correlations were obtained with a 7 .5-ft-diameter 
hingeless rotor model tested in the USAAMRDL— Ames 
wind tunnel.^- The model configuration and test 
conditions covered a wide range of parameters. 

The results Included here are for p * 1.15 and 
advance ratios from 0.0 to 0.6. 

Elastic Blade Bending 

In Fig. 7, experimental values of roll and 
pitch moments due to 0 8 are compared with theo- 
retical results which are calculated neglecting 
induced flow perturbations. Two sets of theory 
are presented. The first theory employs a rigid 
centrally-hinged blade with root spring to model 
the elastic blade bending, and the second theory 
uses a similar model, except that hinge offset is 
allowed. The largest differences between the two 
theories occur near resonant frequencies, i.e., 
a) * 0.15, 1.15. (The primary effect of mode 
shape is aerodynamic, Eq, (8); it causes domi- 
nance at resonance.) A surprising element in 
Figure 7 is that the centrally hinged model gives 
closer agreement with the high frequency response 
than does the hinge offset model. This reversal, 
however, is not a consistent trend in the data 
and may be somewhat coincidental. 



Figure 7. Comparison of experimental data with 
rigid blade approximations without 
induced flow, p * 1.15, y * 4.25, 

B ® 0.97, ep C - 0.25, u - 0.60. 


Similar frequency response comparisons have 
been made when the blade is modeled by one or two 
of the rotating elastic mode shapes. When 
y < 0.8 and w is at least once-per-r evolution 
below the second flap frequency, the one- and two- 
mode calculations are within a few percent of the 
hinge-offset results. At higher advance ratios 
and frequencies, the effects of second-mode bend- 
ing can become significant; but in the range of 
operating conditions considered here, a single ro- 
tating mode is sufficient to model the blade. 


Three major types of discrepancies between 
theory and experiment which are found in Figure 7 


cannot be explained in terras of flapping mode shape 
effects. The first is the difference encountered 
at frequencies near one and two per revolution. 

This difference may be explained by the fact that 
the lead-lag frequency of this configuration is 
near two per revolution, causing resonance at these 
frequencies. The second discrepancy is the irregu- 
larity in the pitch response at w « 0.6. Here, a 
natural frequency of the rotor support stand is 
being excited and contaminates the data.^- The 
third discrepancy is found at u> < 0.6, and will be 
shown to result from unsteady inflow perturbations. 

Effect of Induced Flow In Hover 


The low-frequency hover data provide some in- 
sight into the effects of unsteady induced flow. 

In Figure 8, rotor roll and pitch moments versus 
S g are presented. The experimental results are 
for $ * 4°, v * 0.03. The theoretical results are 
calculated using the actual blade rotating mode 
shape as a generalized coordinate and using three 
different representations of the induced flow. The 
first representation is the elementary model, which 
completely neglects induced flow perturbations. 

The second representation is quasi-steady momentum 
theory, which neglects the apparent inertia 
(Kj * 0), assuming that nonuniform induced flow 
perturbations instantaneously follow the blade dy- 
namics. The third representation is unsteady mo- 
mentum theory, which gives a time lag on the in- 
duced flow perturbations. (The empirical model is 
not applicable in hover.) 

A comparison of theory and experiment reveals 
that the elementary theory is unsatisfactory below 
u » 0.6, failing to reproduce even the qualitative 
character of the data. On the other hand, the 
theories which Include induced flow perturbations 
account for most of the important features of the 
response. The loss of aerodynamic effectiveness, 
which is a result of induced flow perturbations, 
causes a decrease in the excitation forces and an 
overall decrease in the response. But the loss of 
aerodynamic effectiveness also lowers the blade 
damping, causing a resonant peak effect near the 
blade natural frequency (with p * 1.15, u * 0.15). 

The effect of the unsteady induced flow terms 
is also evidenced in Figure 8. The major contri- 
bution of Rj is the determination of how rapidly 
with u) the aerodynamic effectiveness returns to 
the elementary value. Above u> * 0.6, the theo- 
retical value of Kj gives the proper amplitude 
and phase for the hub moments, while the quasi- 
steady theory (K x • 0) fails to return to the con- 
ventional value and does not agree with the data. 
Below a) - 0.6 the comparison is less clear. In 
the roll-moment phase and amplitude, a Kj less 
than 0.1132 would give better correlation than 
does this theoretical value. In the pitch-moment 
response, however, a smaller Kj would give worse 
correlation than does Kj *» 0.1132. Further work 
would be necessary to determine if this effect is 
due to experimental difficulties (such as recir- 
culation) or to an actual deficiency in the in- 
duced flow model. 
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Rotor response to cyclic pitch in ~ 
hover, p - 1.15, y - 4.25, B - 0.97, 
e pc * 0.25, y = 0, aa * 0.7294, 
v e X « 0.03, momentum theory, single 
rotating mode. 


In Figure 9, rotor roll and pitch moments 
versus a are presented for the same test condi- 
tions as in Figure 8. Data are presented for 
shaft excitations in both roll and pitch, since 
in hover the response to these controls is 
ideally symmetric. A comparison of the two sets 
of data gives an indication of the experimental 
error due to test stand dynamics (and possibly 
recirculation). Although the data are question- 
able for a) > 0.3, the lower frequency data sub- 
stantiate three of the observations made from 
Figure 8. First, the elementary theory is quali- 
tatively inaccurate for amplitude and phase re- 
sponse. Second, a major effect of induced flow 
is a resonant peak effect near a) - 0.15. Third, 
Kj < 0.1132 would give better correlation than 
the theoretical value at low u. Figure 9 also 




Figure 9. Rotor response to hub motions in hover, 
p - 1.15, y - 4.25, B - 0.97, 

•pc * 0.25, y - 0, ca - 0.7294, 
v ■ J ■ 0.03, moment tun theory, single 
rotating mode. 


shows that although induced flow decreases the 
blade damping, it can actually increase the rotor 
pitch/rate damping - -Re[3(C M /aa)/3a] and also 
Increase the rotor pitch/roll coupling 
* -Re[3 (CL/aa)/3a] . The damping and coupling can 
be found by dividing the plotted curves by 
-iu), A - iuxx, which is approximately equivalent to 
taking the slope of the plotted curves with a 90- 
degree shift in phase angle. For this particular 
configuration, the damping and coupling are in- 
creased by induced flow effects, indicating that 
induced flow perturbations can be important in 
coupled rotor /fuselage dynamics. 

Effect of Induced Flow in Forward Flight 


In the next three figures, experimental data 
at high-advance ratio (y« 0.51) and very low lift 
(6 ~ 0.5°) are compared with theory using three 
Induced flow descriptions. The first description 
is an analysis which neglects induced flow pertur- 
bations, the second description is the empirical 
model of Reference 2 with no time lag (quasi- 
steady, Kj - Kjj « 0) , and the third description 
is the empirical model of Reference 2 adapted to 
the unsteady case according to Eq. (34) (with the 



Figure 10. Rotor response to collective pitch in 
forward flight, p - 1.15, y - 4.25, 

B - 0.97, e pc - 0.25, y - 0.51, 
aa - 0.7294, v - X - 0, single 
rotating mode. 


theoretical values of Kj and 1^) . The first 
comparison of theory and experiment is shown in 
Figure 10 for the roll- and pitch-moment response 
due to 0 O . The elementary theory predicts a 
roll moment of 0.017 at u ■ 0 and a near-zero 
crossing (amplitude - 0, phase angle discontinu- 
ous) at u> ■ 0*4. The data, however, displays a 
much lower steady value and completely avoids the 
zero crossing. The unsteady and quasi-steady 
empirical models provide a fairly accurate 
description of this behavior, showing quantita- 
tive agreement with phase and magnitude for 
w < 0.6. For the pitch moment derivative, the 
empirical models predict the qualitative (but not 
the quantitative) aspects of the reduction in 
moment (from the conventional value) due to 
induced flow. 
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Another comparison of theory and experiment 
is shown in Figure 11 for the roll- and pitch- 
moment response due to 0 B . The empirical models 
predict a roll-moment derivative which is less 
than the elementary value, exhibiting a near-zero 
crossing at a * 0.26. This characteristic is 
clearly evident in the magnitude and phase of the 
data, but it does not appear in the theory 



W (i) 



Figure 11. Rotor response to longitudinal cyclic 
pitch in forward flight, p ■ 1.15, 
y = 4.25, B » 0.97, ep C * 0.25, 

\i = 0.51, aa - 0.7294, v * X * 0, 
single rotating mode. 


without induced flow. For the pitch-moment deriva- 
tive, the elementary theory agrees with the data 
only for w > 1.2; the quasi-steady theory shows 
good correlation for 0 < to < 0.6, and the unsteady 
theory gives quantitative correlation at all fre- 
quencies. 


The third comparison is shown in Figure 12 for 
the roll- and pitch-moment response due to 0 C . 

The data show that the roll-moment derivative is 
less than the elementary value at w ■ 0, display- 
ing a resonant peak (near to * 0.15) which is 
greater than the elementary value and which is ac- 
companied by a 10-degree phase shift. The empiri- 
cal models predict the qualitative character of the 
resonant peak and quantitative character of the 
phase shift. The empirical models also correlate 
well with the pitch-moment response, for which the 
experiment shows the derivative to be greater than 
the elementary value for a> < 0.3 and less than 
the elementary value for w > 0.3. 

In general, the empirical inflow models show 
this same degree of correlation at all advance 
ratios considered (y “ 0.27, 0.36, 0.51, 0.60). 
This substantiates one of the qualitative 
conclusions of Figure 6. For moderate advance 
ratios and w < 1.0, an appropriate unsteady or 
quasi-steady induced flow theory is adequate, but 
the theory without induced flow is in consider- 
able error. Of course, Figure 6 only Implies in 
which regions quasi-steady or unsteady terms may 
be significant. It does not imply that any par- 
ticular quasi-steady or unsteady model will be 




Figure 12. Rotor response to lateral cyclic 
pitch in forward flight, p = 1.15, 
y « 4.25, B * 0.97, ep C = 0^25, 
p * 0.51, aa - 0.7294, v « X « 0, 
single rotating mode. 

adequate. For example, in Figure 13, pitch mo- 
ment derivates (as calculated using the theory 
without induced flow, unsteady momentum theory, 
and unsteady empirical theory) are compared with 
the experimental data. The comparison shows that 
unsteady momentum theory can be in qualitative 
disagreement with the data even though empirical 
theory shows good correlation. Even the empiri- 
cal model, however, does not show complete quan- 
titative correlation; and further refinements in 
the induced flow model may be necessary. 


Conclusions 


1. On the basis of an equivalent Lock number 
relation and p « 00 , quasi-steady nonuniform in- 
duced flow perturbations can have a significant 
effect on rotor response throughout the entire 
thrust /advance ratio range; but the time lag of 
the induced flow is only important at low lift and 
low advance ratio. 

2. In hover, unsteady momentum theory with appar- 
ent mass terms from potential flow provides a 
significant Improvement in data correlation over 
the theory without induced flow perturbations; 

but further work is required to refine the induced 
flow model. 

3. In forward flight and near-zero lift, the 
empirical inflow model of Reference 2, whether 
used with the unsteady time-lag effect or with- 
out the time-lag effect (quasi-steady) , corre- 
lates well with most qualitative and some quanti- 
tative aspects of the data, while unsteady momen- 
tum theory and the theory without induced flow 
provide little agreement with the data. 

4. A single rotating mode is sufficient for 
flapping response calculations when y <0.8 and 
when the major excitation frequency is at least 
one e-per-r evolution below the second flapping 
frequency. 
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Figure 13. Effect of induced flow model on low frequency, roll response, p * 1.15, y ■ 4.25, B * 0.97, 

e^ c ** 0.25, p = 0.51, Qa * 0.7294, V * X » 0, - 0.1132, » 0.8488, single rotating mode. 
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